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' It is known that the gravitational collapse of a dust ball results in naked singularity formation 

£\J . from an initial density profile which is physically reasonable. In this paper, we show that explosive 

, , . radiation is emitted during the formation process of the naked singularity. 

; PACS numbers: 04.20.Dw, 04.70.Dy, 98.70.Sa 

■ It is known that the gravitational collapse of an inhomogeneous dust ball results in shell-focusing naked singularity 
formation It has been also shown that the naked singularity formation is possible from the spherical collapse 

y—i • of a perfect fluid with a very soft equation of state |^-^]. Moreover, a kind of runaway collapse in Newtonian gravity 
is similar to the naked singularity formation in general relativity in many respects. These strongly suggest that the 
gravitational collapse will often involve drastic growth of spacetime curvature outside the event horizon. 

From this point of view, a number of researchers have examined emission during the naked singularity formation. 
In classical theory, Nakamura, Shibata and Nakao Q suggested that the forming naked singularity in the collapse of a 
prolate spheroid may be a strong source of gravitational waves. Recently, Iguchi, Nakao and Harada |l(J and Iguchi, 
Harada and Nakao pl] , p^ ] examined the behavior of nonspherical linear perturbations of the spherical dust collapse. 
They reported rather milder instability in |l2]. 
' Hawking p3[ showed that thermal radiation is emitted from the gravitational collapse to a black hole by quantum 
effects. In the formation of a globally naked singularity, the spacetime curvature grows unboundedly and the strongly 
curved region can be seen by a distant observer, unlike in the formation of a black hole. This fact suggests that the 
forming naked singularity may be a strong source of radiation owing to quantum effects. In this context, Ford and 
Parker jl4| calculated the radiation during the formation of a shell-crossing naked singularity and found that the 
luminosity is finite. This will be because the shell-crossing singularity is very weak. It is known that the shell-focusing 
singularity is stronger than the shell-crossing singularity. Hiscock, Williams and Eardley fl5|| showed the diverging 
luminosity during the formation of the shell-focusing singularity in the spherically symmetric, self-similar collapse 
of a null dust. Barve, Singh, Vaz and Witten []l6|| also showed the diverging luminosity during the shell-focusing 
singularity in the self-similar collapse of a dust ball. See also |T7j , |l8| ]. 

The last two examples in which the diverging luminosity is emitted are self-similar collapse. However, the self-similar 
collapse is a particular solution among gravitational collapse solutions in general relativity. Moreover, it is uncertain 
whether or not the central part of the realistic spherical collapse with nonzero pressure tends to be self-similar in 
strong-gravity regime such as shell- focusing singularity formation (cf. || ) . For the self-similar collapse of a dust ball, 
it has been shown that the redshift at the center diverges to infinity and that the curvature strength of the naked 
singularity is very strong. In fact, it is known that this solution does not allow an initial density profile which is a C°° 
function with respect to the local Cartesian coordinates. For C°° case, the features are much different. The redshift 
is finite and the curvature strength is not very strong [fl9|-pl| . Though Einstein equation does not require such strong 
differentiability of initial data, we usually set such initial data in most astrophysical numerical simulations. We should 
also comment that this model of the collapsing dust ball will be valid until perturbations sufficiently grow due to the 
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reported mild instability. In this paper, we examine radiation during the naked singularity formation in the collapse 
of an inhomogeneous dust ball from an initial density profile which is a C°° function. We use the units in which 
G = c = h=l. 

We consider both minimally and conformally coupled massless scalar fields in four dimensional spacetime which is 
spherically symmetric and asymptotically flat. Let u and v be null coordinates such that they are written as u rs T '— R 
and v ~ T + R in the asymptotic region with the quasi-Minkowskian spherical coordinates (T, R, 9, <fr). An outgoing 
null ray u = const arriving on 3 + can be traced back through the geometry becoming an incoming null ray v = const 
originating from 5~ with v. This gives the relation between u and v and we define the function G(u) by v = G(u). 
Here, we assume that geometrical optics approximation is valid, which implies that the trajectory of the null ray gives 
a surface of a constant phase of the scalar field. Then, the luminosity Li m for the minimally coupled scalar field and 
the luminosity L\ m for the conformally coupled scalar field for fixed I and m are given through the point-splitting 
regularization as 

T 1 f9lY _ 1 f G "\ f I (G"V m 

lm ~ 48tt \G'J 24tt \G' J ' lm ~ 48tt \G' J ' 1 ' 

This implies that the luminosity depends on how the scalar field couples with gravity. However, if G" /G'\ u — a = 
G" /G'\ u= b holds, the amounts of the radiated energy during a < u < b for the both fields are the same. We should 
note that the geometrical optics approximation is only valid for smaller Hereafter we omit the suffix I and m. It is 
noted that these results are free of ambiguity coming from local curvature because the regularization is done in flat 
spacetime. 

The spherically symmetric collapse of a dust fluid is exactly solved (2|j2^]. The solution is called the Lemaitre- 
Tolman-Bondi (LTB) solution. For simplicity, we assume the marginally bound collapse. The metric is given by 

ds 2 = -dt 2 + R 2 r (t, r)dr 2 + R 2 {t, r)dfl 2 (2) 

in the synchronous comoving coordinates with dH, 2 = d9 2 + sin 2 9dcj} 2 . The energy density is given by 

, F'(r) 
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where F(r) is equal to twice the Misner-Sharp mass. We rescale the radial coordinate r as R(Q, r) = r. Then, R is 
given as 

/ / — \ 2/3 
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The singularity occurs at the time t s (r) = (2/3)y/r 3 /F. We denote the time of occurrence of singularity at the center 
as to = t s (0). From Eq. (||), if we require that an initial density profile at t = is a C°° function with respect to the 
local Cartesian coordinates, F(r) is expanded around r = as 

F(r) = F 3 r 3 + F 5 r 5 + F 7 r 7 + ■■■, (5) 

where we assume that F$ is positive. If we assume the marginally bound collapse and the mass function as Eq. (|^), 
then the central shell- focusing singularity is naked if and only if F 5 is negative p]-^,p4]-p6|| . At an arbitrary radius 
r = r s f, the LTB spacetime can be matched with the Schwarzschild spacetime 



ds 2 = - [1- — ) dT 2 + [I- — ) dR 2 + R 2 dfl 2 , (6) 
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with M = F(r s f)/2. The relation between T and t is given at r = r s f as 



T = t = - 2V2MR + 2M In V I_ V ; (7) 
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The null coordinates u and v in the Schwarzschild spacetime are defined as 

u = T-R*, v = T + R* (8) 
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with R* = R + 2M\n[(R/2M) - 1]. 

We can determine the function G by solving the trajectories of outgoing and ingoing null rays in the dust and 
calculating u and v through Eq. (||) at the time when the outgoing and ingoing null rays reach the surface boundary. 
The trajectories of null rays in the dust are given by the following ordinary differential equation 

dt „ 

— = ±R, r , (9) 
ar 

where the upper and lower signs denote outgoing and ingoing null rays, respectively. We have numerically solved 
Eq. (S) by the Runge-Kutta method of the fourth order. We have carried out the quadruple precision calculation for 
retaining accuracy. We have chosen the mass function as F(r) — F^r 3, + F$r 5 . We hnd that the central singularity 
is globally naked for very small r s f if we fix the value of F3 and F5. Although we have calculated several models, 
we only display the numerical results for the model with F3 = 1, F5 = —2 and r s t — 0.02 in an arbitrary unit 
because the features are the same if the singularity is globally naked. The total gravitational mass M is given by 
M = 3.9968 x 10~ 6 for this model. 

See Fig. [|. We define uo a s the retarded time of the earliest light ray which originates from the singularity. In 
Fig. [l](a), the first derivative G'(u) is plotted. This shows that G'(u) does not diverge but converge to some positive 
value A with < A < 1. In Fig. [j](b), it is found that the second derivative G"{u) does diverge as in no. This figure 
shows that the behaviors of growth of G"(u) are different each other for 10~ 4 < uq — u and for < uo — u ^, 10~ 4 . For 
10~ 4 < u — u , the dependence on u is written as G" oc — (— u)~ 2 , while, for < u — u < 10~ 4 , the dependence is 
written as G" oc — (uq — u) -1 / 2 . Here we determine a dimensionful constant of proportion by physical consideration. 
We assume that the early time behavior is due to the collapse of the dust as a whole while the late time behavior 
is due to the growth of the central curvature. This assumption will be justified later. First we should note that the 
coefficient must be written using the initial data because it must not depend on time. Next we should note that we 
can regard any t — const < to hypersurface as an initial hypersurface. Therefore, the coefficient must be independent 
of the choice of an initial slice. For the early time behavior, the only possible quantity is the gravitational mass M 
of the dust cloud. For the late time behavior, the only possible quantity is lu s = igAo = (3/2) 7 F 3 13 / 2 (— Fs) -3 , where 

to is given as to = (2/3)F^ 1 ^ 2 and Iq denotes the scale of inhomogeneity defined as Iq = (— .F5/.F3) -1 / 2 . We call uj s a 
singularity frequency. The lu s is independent of the choice of an initial slice because the mass function F(r) is written 
in terms of R as 

^0 ~ t\ 7->3 , {to ~ 



F(r) = F 3 l-^j— \ r' + f 5 \^— j R J + ---, (10) 

around the center. Thus, we can write G" ~ ~f e M(—u)~ 2 for — u M, and G" ~ — fiAioV 2 '(uq — u)^ 1 / 2 for 
< uq — u <C uj 1 , where f e and fi are dimensionless positive constants of order unity. The turning point from the 
early time behavior to the late time behavior is roughly estimated as uq — u w (Mujs) 2 ^ 3 ^ 1 . These estimates show 
a good agreement with the numerical results. 

Let us consider the luminosity which is calculated by Eq. (|l]). The numerical results are displayed in Fig. ||. We 
can also write an analytic expression for the luminosity using Eq. ([!]). In Fig. 0(a), we can find G' ~ A for the late 
time behavior. Then, the luminosity for the late time behavior is obtained as 

L « ^/^K - «r 3/2 , L « ^/j 2 w.(uo - u)" 1 . (11) 

Therefore the luminosity diverges to positive infinity for the both fields as u — > uq. The radiated energy is obtained 
by integrating the luminosity with u. The radiated energy for the late time behavior is estimated as 

E^^T(u Q -u)-^, £~-L/-Wn {M " s?/ \ . (12) 

247T 4»7T OJ s {Uo — U) 

Therefore, the amounts of the total radiated energy diverge to positive infinity for the both fields as u — > uq. In 
realistic situations, we may assume that the naked singularity formation is prevented by some mechanism and that 
the quantum particle creation is ceased at the time Uq — u « At, which implies that G" /G' vanishes for u — u < At. 
Then, the second term in the expression of the luminosity L gives no contribution to the total radiated energy. 
Therefore, the amounts of the total energy for the both fields are the same, i.e., 

E = E^^ff.M^f\ (13) 
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From the numerical results and physical discussion above, we obtain the following formula for the late time behavior 



G(u) w A(u - u Q ) - -Afiu l J 2 (u - u) 3/2 + const. (14) 
o 

For comparison and for a test of our numerical code, we have also calculated the function G for the Oppenheimer- 
Snyder (OS) collapse to a black hole which is given by F(r) = F^r 3 . We only display the numerical results for the 
model with F3 = 1 and r s f = 0.02 in an arbitrary unit. M is given by M = 4 x 10~ 6 for this model. The results are 
plotted also in Figs, [j] and |[ Because we cannot define Uq for the OS spacetime, we plot the numerical results by 
setting uq — 0. For u > 0, the numerical results show the well-known behavior 

G(u) fa -const • cxp ^j^J + v h , (15) 

LR ^ 768^ (16) 

where the ingoing null ray with Vh is reflected to the outgoing null ray which is on the event horizon. The numerical 
results of L and L for the OS collapse show a good agreement with Eq. (|l(|). In Figs. |l| and |], we can find that the 
assumption concerning the early time and late time behaviors is justified. 

If a back reaction of quantum effects would not become important until a considerable fraction of the total energy 
of the system is radiated away, the emitted energy could amount to E ~ 1O 54 (M/M )erg. Thus, the naked singularity 
explosion would be a new candidate for a source of ultra high energy cosmic rays. It may also be a candidate for the 
central engine of a gamma ray burst. In order to study such possibilities, we need to transform the obtained intrinsic 
quantities to observed quantities by taking possible reactions of created energetic particles into consideration. It is 
noted that, since we have obtained the formula for the late time behavior of the function of G(u), we can determine 
the spectrum of the radiation and thereby estimate the validity of the geometrical optics approximation. We are now 
obtaining positive evidences for the consistency |27j . 
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FIG. 2. Luminosity for minimally and conformally coupled scalar fields in the LTB spacetime and in the OS spacetime 
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